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Abstract 



We present a novel algorithm for deciding whether a given planar curve is an image of a 
given spatial curve, obtained by a central or a parallel projection with unknown parameters. 
A straightforward approach to this problem consists of setting up a system of conditions on 
the projection parameters and then checking whether or not this system has a solution. The 
main advantage of the algorithm presented here, in comparison to algorithms based on the 
straightforward approach, lies in a significant reduction of a number of real parameters that 
need to be eliminated in order to establish existence or non-existence of a projection that maps 
a given spatial curve to a given planar curve. Our algorithm is based on projection criteria 
that reduce the projection problem to a certain modification of the equivalence problem of 
planar curves under affine and projective transformations. The latter problem is then solved by 
differential signature construction based on Cartan's moving frame method. A similar approach 
can be used to decide whether a given finite set of ordered points on a plane is an image of a 
given finite set of ordered points in M.^. The motivation comes from the problem of establishing 
a correspondence between an object and an image, taken by a camera with unknown position 
and parameters. 

Keywords: Central and parallel projections, finite and affine cameras, curves, separating differ- 
ential invariants, projective and affine transformations, geometric invariants, signatures, machine 
vision 
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1 Introduction 



A central projection from to M? models a simple pinhole camera pictured in Figure 1 . A generic 
central projection is described by a linear fractional transformation ([T]): 




PU Zl + Pl2 Z2 + Pl3 Z3 + Pli 
P31 Zl + P32 Z2 + P33 Z3 + P34 



P2I Zl + P22 Z2 + P23 Z3 + P24 



(1) 



^ P31 Zl + P32 Z2 + P33 Z3 + P34 

Figure 1. 

where (2:1,^2,-23) denote coordinates in M^,(x,y) denote coordinates in and pij, i = 1 . . .3, j = 
1 ... 4, are real parameters of the projection, such that the left 3x3 submatrix of 3 x 4 matrix 
P = {pij) has a non zero determinant. The parameters represent a freedom to choose the center of 
a projection, a position of image plane and, in general, non-orthogonal coordinate system on the 
image planej^ In the case when the distance between a camera and an object is significantly greater 
than the object depth a parallel projection provides a good camera model. A parallel projection 
has 8 parameters and can be described by a 3 x 4 matrix of rank 3, whose last row is (0, 0, 0, 1). We 
review various camera models and related geometry in Section [2] (see also ||14|). In most general 
terms, the projection problem considered here can be formulated as follows: 

Problem 1. Given a subset Z o/M'^ and a subset X of^ , does there exist a projection P : — )• 
such that X = P{Z) ? 

A straightforward approach to this problem consists of setting up a system of conditions on the 
projection parameters and then checking whether or not this system has a solution. In the case when 
Z and X are finite lists of points, a solution based on the straightforward approach can be found 
in [14]. For curves and surfaces under central projections, this approach is taken in [9j. However, 
internal parameters of the camera are considered to be known in that paper and, therefore, there 
are only 6 in that study vs. 12 free parameters considered here. The method presented in (9] also 
uses an additional assumption that a planar curve C has at least two points, whose tangent 
lines coincide. An alternative approach to the problem in the case when Z and X are finite lists 
of points under parallel projections was presented in jl||2]. In these articles, the authors establish 
polynomial relationships that have to be satisfied by coordinates of the points in the sets Z and X 
in order for a projection to exists. Our approach to the projection problem for curves is closer in 
spirit to (l]|2], as we also establish necessary and sufficient conditions that curves Z and X must 
satisfy in order for a projection to exist. However, unlike [l|[2], we exploit the relationship between 
the projection problem and equivalence problem under group-actions. We will show below that, in 
comparison with the straightforward approach, our solution leads to a significant reduction of the 
number of parameters that have to be eliminated in order to solve Problem [T] for curves. 

Throughout the paper, we assume that Z and X are algebraic curves and Z is not a straight line 
(and, therefore, its image under any projection is a one-dimensional constructible set). We will also 
relax the projection condition to X = P{Z), where bar throughout the paper denotes the algebraic 



^It is clear from ([l]) that multiplication of P by a non-zero constant does not change the projection map, therefore, 
we can identify P with a point of the projective space PR^^, rather then a point in R^^. However, since we do not 
know which of the parameters are non-zero, in computations we have to keep all 12 parameters. 
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closure of a set. In this section, as well as in SectionjSj where algorithms and examples are presented, 
we make an additional assumption that Z and X have rational parameterization: Z = {T{s)\s G M} 
and X = {7(t)|t € M}, where T: M and 7: M are rational maps. This assumption is not 

essential. Algorithms for non-rational algebraic curves are placed in Appendix 8.5 Problem [T] for 
central projections of rational algebraic curves can be reformulated as the following real quantifier 
elimination problem: 

Reformulation 1. ("straightforward approach) Given rational mapsT: R — )• andj: R 
R^, determine the truth of the statement: 

3PGR3X4 det(pij)&3V Vs € R 3t G R -f{t) = P{T{s)). 

Real quantifier elimination problems are algorithmically solvable ||22j. A survey of subsequent 
developments in this area can be found, for instance, in 16 and [g]. High computational complexity 
of these problems makes a reduction in the number of parameters to be desirable. 

The projection criteria developed in this paper reduces the projection problem to the problem 
of deciding whether the given planar curve X is equivalent to a curve in a certain family of planar 
curves under an action of the projective group in the case of central projections, and under an 
action of the affine group in the case of parallel projections. The family of curves depends on 3 
parameters in the case of central projections, and on 2 parameters in the case of parallel projections. 

These group-equivalence problems can be solved by an adaptation of differential signature con- 
struction developed in [s] to solve local equivalence problems for smooth curves. In Section |4| 
we give an algebraic formulation of signature construction and show that it leads to a solution of 
global equivalence problems for algebraic curves. For this purpose we introduce a notion of differen- 
tially separating set of invariants. Following this method for the case of central projections, when 
Z and X are rational algebraic curves, we define two rational signature maps : R R^ and 
Sz : R'* — )• R^ . Construction of these signature maps requires only differentiation and arithmetic 
operations and is computationally trivial. Problem [l] for central projections of rational algebraic 
curves becomes equivalent to: 

Reformulation 2. ("signature approach) Given two rational maps Sx ■ ^ ^ and : R"^ — ^ 

R^, determine the truth of the statement: 3c G R'^ Vs G R 3t G R Sz{c, s) = Sx{t). 

We note that Reformulation [T] and Reformulation [2] have similar structure, but the former 
requires elimination of 14 parameters, while the latter requires elimination of only 5 parameters. 
(The case of parallel projection is treated in the similar manner and leads to the reduction of the 
number of real parameters that need to be eliminated from 10 to 4.) 

Another advantage of the approach presented here is its universality: essentially the same 
method can be adapted to various types of the projections and various types of objects, both 
continuous and discrete. In Section [6| we discuss how this method leads to an alternative solution 
to the projection problem for finite lists of points (see [4j for more details). Our method can be also 
easily adapted to solve projection problems from R*^ to a hyperplane. The problem and the solution 
remains valid if we replace R with C, and in fact, as with most problems in algebraic geometry, the 
implementation of the algorithms is easier over the complex numbers (see also Remark 21 below). 
The existence of a complex projection provides a necessary but not a sufficient condition for the 
existence of a real projection. 

Although the relation between projections and group actions is known, our literature search 
did not yield algorithms that exploit this relationship to solve the projection problem for curves 
in the generic setting of cameras with unknown internal and external parameters. The goal of 
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the paper is to introduce such algorithms. Although the development of efficient implementations 
of these algorithms and their complexity study lie outside of the scope of this paper, we made 
a preliminary implementations in Maple of projection algorithms over complex numbers. We 
implemented both an algorithm based on signature construction presented here and an algorithm 
based on the straightforward approach and timed their performance. The code and the experiments 
are posted at www.niath.ncsu.edu/-iakogan/symbolic/projections.html, 

In order to become practically useful in real-life applications, the algorithmic solution presented 
here, would have to be adapted to curves given by finite sampling of points. Some directions of 
such adaptation are indicated in Section [7] of the paper. 

2 Projections and cameras 

We will embed M" into the projective space PR" and use homogeneous coordinates on PM*^ to 
express the map ([T]) by matrix multiplication. 

Notation 1. Square brackets around matrices (and, in particular, vectors) will be used to denote 
an equivalence class with respect to multiplication of a matrix by a nonzero scalar. Multiplication 
of equivalence classes of matrices A and B of appropriate sizes is well-defined by [A] [B] := [AB]. 

With this notation, a point (x, y) G corresponds to a point [x, y, 1] = [Ax, Ay, A] G PM^ for 
all A 7^ 0, and a point {zi,Z2,zs) S corresponds to [zi, Z2, zs,l] G PM^. We will refer to the 
points in PM" whose last homogeneous coordinate is zero as points at infinity. In homogeneous 
coordinates projection ([l]) is a map [P] : PM'^ — t- PM^ given by 

[x,y,l]^ = [P][zu Z2, zs,lf, (2) 

where P is 3 x 4 matrix of rank 3 and superscript T denotes transposition. Matrix P has a 1- 
dimensional kernel. Therefore, there exists a point [zj* , , , ^4] G PM'^ whose image under the 
projection is undefined (recall that [0,0,0] is not a point in PM^). Geometrically, this point is the 
center of the projection. 

In computer science literature (e.g. [14]), a camera is called finite if its center is not at inffiiity. A 
finite camera is modeled by a matrix P, whose left 3x3 submatrix is non-singular. Geometrically, 
finite cameras correspond to central projections from to a plane. On the contrary, an infinite 
camera has its center at an infinite point of PM^. An infinite camera is modeled by a matrix P 
whose left 3x3 submatrix is singular. An infinite camera is called affine if the preimage of the 
line at infinity in PM^ is the plane at infinity in PM'^. An affine camera is modeled by a matrix P 
whose last row is (0, 0, 0, 1). In this case map ([l]) becomes x = pn zi + pu Z2 + P13 z^ + pi^ , y = 
P21 z\ + P22 Z2 +P23 -23 + P2A- Geometrically, affine cameras correspond to parallel projections from 

to a planej^ Eight degrees of freedom reflect a choice of the direction of a projection, a position 
of the image plane and a choice of linear system of coordinates on the image plane. An image 
plane may be assumed to be perpendicular to the direction of the projection, since other choices 
are absorbed in the freedom to choose, in general, non-orthogonal coordinate system on the image 
plane. 

Definition 2. The set of equivalence classes [P], where P = (j'ij)j=\'"4 is a 3 x 4 matrix whose 
left 3x3 submatrix is non-singular, is called the set of central projections and is denoted CV. 

^Parallel projections are also called generalized weak perspective projections lll sl. 
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The set of equivalence classes [P], where P = {pijYj^i''^ has rank 3 and its last row is (0, 0, 0, A), 
A 7^ 0, is called the set of parallel projections and is denoted VV. 

Equation ([T]) determines a central projection when [P] € CV and it determines a parallel pro- 
jection when [P] £ VV. Sets CV and VV are disjoint. Projections that are not included in these 
two classes correspond to infinite, non-affine cameras. These are not frequently used in computer 
vision and are not considered in this paper. 



3 Projection criteria for curves 

Recall that for every algebraic curve X C M" there exists a unique projective algebraic curve 
[X] C PM" such that [X] is the smallest projective variety containing X (see ||,iT] ) . 

Definition 3. We say that a curve 2 C M'^ projects to A' C if there exists a 3 x 4 matrix P of 
rank 3 such that [X] = {[P][z] | z G 2, Pz / 0}. We then write X = P{Z) or [X] = [P][Z]. 

Definition 4. The projective group VQC{n + l) is a quotient of the general linear group QL[n + l), 
consisting of (n + l) x (n + 1) non-singular matrices, by a 1-dimensional ahelian subgroup XI, where 
A 7^ G M and I is the identity matrix. Elements ofVQC{n + l) are equivalence classes [B] = [\B], 
where A 7^ and B G QL{n + 1). 

The affine group A{n) is a subgroup of VQC{n + 1) whose elements [B] have a representative 
B G QL{n + 1) with the last row equal to (0, ... ,0, 1). 

The special affine group SA{n) is a subgroup of A{n) whose elements [B] have a representative 
B G QL{n + 1) with determinant 1 and the last row equal to (0, . . . , 0, 1). 

In homogeneous coordinates, the standard action of the projective group VQC{n + 1) on PR"" 
is defined by multiplication: 

[Zi, ...,Zn, Zq]^ [B][zi, . . . , Zn, Zq]^ . (3) 

The action ([s]) induces linear-fractional action of VQC{n + 1) on M". In particular, for n = 2, 
[B\ G VQC{?>) we have 

(xy\^( ^11 x + bi2y + bi3 hi X + b22y + b23 \ .^x 

\b31x + b32y + &33 ' b3ix + 632 y + bssj ■ 

The restriction of ([s]) to A{n) induces an action on consisting of compositions of linear trans- 
formations and translations. In particular, for n = 2 and [B] G ^(2) represented by a matrix B 
whose last row is (0, 0, 0, 1) 

(6iix + 6i2y + 6i3, ^21 2; + 622 y + ^23) ■ (5) 

Definition 5. We say that two curves Xi C M" and X2 C M" are VQC{n + l)-equivalent if there 
exists [A] G VgC{n + 1), such that [X2] = {[A][p] \ [p] G [^"1]}. We then write X2 = A{Xi) or 
[X2] = [A][Xi]. If [A] G G, where G is a subgroup ofVQC{n + 1), we say that Xi and X2 are 
G -equivalent. 

Before stating the projection criteria, we make the following simple, but important observations. 
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Proposition 6. (i) If Z C M.^ projects to A' C by a parallel projection, then any curve 
that is A{3)- equivalent to Z projects to any curve that is A{2)- equivalent to X by a parallel 
projection. In other words, parallel projections are defined on affine equivalence classes of 
curves. 

(a) If Z <Z ^ projects to X <Z ^ by a central projection then any curve in M'^ that is .4(3)- 
equivalent to Z projects to any curve on that is VQ C{?>)- equivalent to X by a central 
projection. 

Proof, (i) Assume that there exists a parallel projection \P] G VV such that [X] = [P][Z]. Then 
for all {A, B) G ^(2) x ^(3) we have [A] [X] = [A] [P] [p-^] [B] [Z]. Since [A] [P] [p-'^] G VV , curve 
B[Z) projects to A{X). (ii) is proved similarly. □ 

Remark 7. It is known that if Xi and X2 are images of a curve Z under two central projections 
with the same center, then Xi and X2 are VGC{3)- equivalent, but if the centers of the projections 



are not the same this is no longer true (see Example 28). Similarly, images of Z under different 



parallel projections may not be A{2)- equivalent (see Example 31). 



Theorem 8. ^central projection criteria^ A curve Z C projects to a curve X C R'^ by 
a central projection if and only if there exist ci, C2, C3 G M such that X is VGC{3)- equivalent to 



-'Cl,C2,C3 



Zi + Cl Z2 + C2 

Z3 + C3' Z3 + C3 



(^1, 2^2, 2:3) e z} c 



(6) 



Proof. 



Assume there exists a central projection [P] such that [X] = [P] [Z]. Then P is a 



3x4 matrix P = {pij)^jZ}i"\ whose left 3x3 submatrix is non-singular. Therefore there exist 
Cl, C2, C3 G M such that = cip*i -|- C2P*2 + £3^*3, where p^:j denotes the j-th column of the 
matrix P. We observe that 

(7) 



[A][P^][B] = [P] , 



where A := (p. 



VJj=1...3 



^ is 3 X 3 submatrix of P, 




and B 



( 1 



Vo 



C2 

C3 
1 / 



(8) 



{Z\ 01,2:2 + 02,23 + C3] 



Note that \A\ belongs to VgC{3). Since 

[P^][B][ZI,Z2,Z3, if 

then [X] = [A][Z,, 

,02,03] J where ^ci,c2,c3 is defined by ^ 
(<;=) To prove the converse direction we assume that there exists [A] G VGC{3) and ci, C2, C3 G 
such that [X] = [A][Zci^c2,c3], where [Zci,c2,c3] is defined by ([6). A direct computation shows that 
Z is projected to X by the central projection [P] = [A] [P^] 



B], where B and [P^] are given by 



(|8J). □ 

Theorem 9. (^parallel projection criteria.) A curve Z C R^ projects to a curve X C R^ 
by an parallel projection if and only if there exist ci,C2 G R and an ordered triplet {i,j,k) G 
{(1, 2, 3), (1, 3, 2), (2, 3, 1)} such that X is A{2) -equivalent to 



C\,C2 



{( 



Zi Cl Zfc, 



Zj+C2Z-k) (21,22,23) G C 



(9) 
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Proof. (=^)Assume there exists a parallel projection [P] such that [X] = [P] [Z]. Then [P] can be 
represented by a matrix: 

/ Pu Pl2 Pl3 PlA \ 

-P = P2l P22 P23 P2A (10) 

\ 1 / 

of rank 3. Therefore there exist 1 < i < ? < 3 such that the rank of the submatrix ( ^^-^ ) 

is 2. Then for 1 < A; < 3, such that k ^ i and k ^ j, there exist ci, C2 G M, such that 



Pik 

P2k 



Pit Pij Pm 



ci ( I+C2 ( ^^■^ )• We define j4 := I poi P2i P2A I and define S to be the matrix whose 
\P2i J \ P2j J \ 1 

columns are vectors 6*i := (1,0,0,0)"'", 6*j := (0,1,0,0)^, 6*^ := (ci, C2, 1, 0)""", b^4 = (0,0,0,1)"'". 

/ 1 \ 

We observe that = [P], where po := 1 . Since [P^][B][Z] = [zHt], then 

V 1 / 

[X] = [A][Zc^^c2]- Observe that [A] G ^(2) and the direct statement is proved. 

(<^) To prove the converse direction we assume that there exist [A] £ .4(2), two real numbers ci 
and C2, and a triplet of indices {i, j, k) £ {(1,_2, 3), (1,3,2), (2,3, 1)}, such that [X] = [A][Zc^^'c2], 



where the planar curve Zc^'^c2{s) is given by ([9|). Let Pp and i? be a matrix defined in the first 
part of the proof. A direct computation shows that Z is projected to X by the parallel projection 
[P\ = [A][P^][B\. □ 

The families of curves Zc^^c2 given by have a large overlap. The following corollary eliminates 
this redundancy and, therefore, is useful for practical computations. 

Corollary 10. (reduced parallel projection criteria) A curve Z C projects to X cM? 
by a parallel projection if and only if there exist oi , 02, 6 G M such that the curve X is A{2) -equivalent 
to one of the following planar curves: 



2ai,a2 = 1(21+01^^3,^2 + 02^3) (zi, Z2, Z^) £ Z^ C R'^ , (11) 



2b = i {zi + hz2, z-i 
Z ■ 



{zi,Z2,Z3)eZ^ CR^ (12) 



{(Z2, Z3) |(zi,Z2,Z3)G2}cM2. (13) 



for any ci, C2 G M the set Zl'^^'c^ = "{ (-^i + ci -Zfc, zj + C2 z^ 



Proof. We first prove that for any permutation {i,j,k) of numbers (1,2,3), such that i < j, and 

zi,Z2,Z3) G Z^ is ^(2)-equivalent to 
one of the sets listed in (llll) - (|13l). 

~ 1 2 3 

Obviously, 2^ci,c2 = ^ai,a2 with ai = ci and 02 = C2. 

For Zl^^i, if C2 / then (I ) = f 7 I'' 1 and so Z^^^^ is ^(2)- 

V i 7 V ^3 + C2Z2 J \Z2 + ^^Z3 J 

equivalent to 2^ai,a2 with ai = — ^ and 02 = — . Otherwise, if C2 = 0, the 2^ci,c2 = with 
b = ci. 
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""2 3 1 ~2 3 1 ~ 1 

Similarly for 2^ci,c2 ) if C2 7^ then 2^ci,c2 is ^(2)-equivalent to 2^ai,a2 with ai = — and 02 



Otherwise, if C2 = 
with h = otherwise ci 



then ^ci^c2 = (-22 + ci^i, Z3). If ci / then -^c^^ 
and ^ci% = ^• 



.£1 

C2 ' 



is ^(2)-equivalent to Zi, 



We can reverse the argument and show that any curve given by ( 11 )-(|l3[) is ^(2)-equivalent to 



a curve from family ([9]). Then the reduced criteria follows from Theoremjoj 



□ 



4 Group-equivalence problem 

Theorems [S] and [9] reduce the projection problem to the problem of establishing group-action equiv- 
alence between a given curve and a curve from a certain family. A variety of methods exist to solve 
group-equivalence problem for curves. We base our algorithm on the differential signature con- 
struction described in [s]. The differential signature proposed there solves the local equivalence 
problem for smooth curves. We adapt this construction to the algebraic setting and prove that 
the differential signature gives a solution of the global equivalence problem in the case of algebraic 
curves. For this purpose, we introduce a notion of a differentially separating set of invariants. We 
discuss the possibility of using some other methods for solving the equivalence problem in Section [7| 

4.1 Differential invariants for planar curves 

A rational action of an algebraic group G on can be prolonged to an action on the n-th jet space 
jn _ coordinates (x, y, y^^\ . . . , y*-"^) as follows]^ For a fixed g G G, let (x, y) = g - {x, y). 

Then x, y are rational functions of (x, y) and 

g ■ (x, y, . . . , y^) = (x, y, y^^\ where (14) 

i-.H^^y)y"' ii^i-^y)] 



In (|14|), ^ is the total derivative, applied under assumption that y is function of x\ \ We note that 



the natural projection vr^: J"- ^ J , k < n is equivariant with respect to action (^14|. 

Definition 11. A function on J" is called a differential function. The order of a differential 
function is the maximum value of k such that the function explicitly depends on the variable y^^^ . 



A differential function which is invariant under action is called a differential invariant. 



Remark 12. Due to equivariant property of the projection vr^: J" — t- , k < n, a differential 
invariant of order k on can be viewed as a differential invariant on J" for all n> k. 

Definition 13. Let G act on M^. A set I of rational invariants is separating on a subset W C 
if W is contained in the domain of definition of each I £ I and \/wi, W2 € W 

I{wi) = I{w2), V/ G X <^=^ 3g £ G such that wi = g ■ W2- 



^Uere y = y'^°^ and J° = 

* We note the duality of our view of variables y'*'. On one hand, they are viewed as independent coordinate 
functions on J". On the other hand, the operator ^ is applied under assumption that y is a function of x and, 
therefore, yC^) is also viewed as the fe-th derivative of y with respect to x. (The same duality of view appears in 
calculus of variations.) 
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Definition 14. (^differentially separating set) Let r- dimensional algebraic group G act on 
M^. Let K and T he rational differential invariants of orders r — 1 and r, respectively. The set I = 
{K,T} is called differentially separating if K separates orbits on Zariski open subset W^'~^ C J^'^^ 
andl = {K,T} separates orbits on a Zariski open subset ofW' C J*". 

4.2 Jets of curves and signatures 

Let A' C be defined by an implicit equation F{x, y) = 0. Then the derivatives of y with respect 
to X are rational functions on X , computed using implicit differentiation. Namely, for p G Af: 

(2), , _ -F.Ap) Fyjp) + 2 F^yjp) F.(p) F,(p) - Fyyjp) F|(p) 



P) 



Definition 15. The n-th jet of a curve X cM? is a rational map j^: X ^ J", where for p G A" 

jj(p)=(x(p),y(p),yj)(p),...,yi")(p)). 



From ([T4| it follows that fg-x (5 • P) = 5 • \fx (p)] 
If X has a rational parameterization 7(i) = y{t)), then 



#(7(t)) = i and #(7(t)) ^ 



(fc-i) 



2,... 



(16) 



(17) 



where ' here and below denotes the derivative with respect to the parameter. In this case, we view 
as 3j rational function M — )• J". 



Definition 16. A restriction of a differential function J" 
of <I> with the n-th jet of curve, i. e. = <I> o j J : A" — )• M. 



to a curve X is a composition 



Definition 17. Let X = {i^T, T} he a differentially separating set of invariants for the G-action 
(see Definition 14)- Then a point p £ X is called I -regular if: 

(1) p is a non- singular point of X ; 

(2) j^-i(p) G W'~^ andj^M ^ i 



(3) 



dK 



7^ and 



dT 



i^(p) 



An algebraic curve X C M.'^ is called non-exceptional with respect to I if all but a finite number of 
its points are X-regular. 

Remark 18. If X is a non- exceptional curve with respect to X, then K\x and T\x are rational 
functions on X . The X-regular points of X are in the domain of the definition of these functions. 
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Definition 19. (^signatureJ Let T = {iT, T} he differentially separating set of invariants with 
respect to G-action and X be a non- exceptional curve with respect to I. The signature Sx is the 
image of the rational map S\x : X ^ M? defined by Sx{p) = {K\x{p) , T\x{p))- 



Remark 20. Let I = {K,T} be a set of differentially separating invariants (see Definition I4) 



Let X C be a non X-exceptional curve defined by an implicit equation F{x, y) = 0. Then 

K\xix,y) = kl^(^'y) ' ~ ll\x'y) ' w^^'^^ ki, k2 and ti,t2 are pairs of polynomials with no non- 
constant common factors modulo F. Consider an ideal 

X := {F, k2 >c — ki, t2T — ti, k2t2cr — 1) C M[x, r, x, y, a]. (18) 

The algebraic closure Sx of signature setSx is the variety of the elimination ideal X = XnM[x, r]. 

Remark 21. We note that dvcuSx = if and only if Kx and Tx are constant functions on X 
and dim5;t = 1 otherwise. Ln the latter case Sx is an algebraic planar curve with a single defining 
equation Sx{x,t) = 0. The equality of signatures for two curves, Sxi = Sx2, implies Sxx{x,t) is 
equal up to a constant multiple to S'^'al^', t). The converse, is not true over M, because it is not an 
algebraically closed field, but is true overC (see [7]j. This, in particular, makes an implementation 
of projection algorithms over real numbers more challenging. 

Theorem 22. (group-equivalence criterion) Assume that algebraic curves Xi and X2 are 
non- exceptional with respect to differentially separating invariants X = (K, T) under G-action. 
Then Xi and X2 are G -equivalent if and only if their signatures are equal: 

—G ^2 ^^=^ Sxi = 8x2 ■ 

Direction =^ follows immediately from the definition of invariants. Direction is proved in 



Appendix 8.2 



Remark 23. A definition of signature for sufficiently smooth non-algebraic curves appears, for 
instance, in [s] and is analogous to our Definitional^ (rationality of functions is no longer required). 
In the smooth case, the direction =^ of Theorem^2^is still valid, but the direction is not. For 
example, curves y = cos{x),x E [0,27r], and y = sin(x),x G [0, 27r], have the same Euclidean 
signatures, but are not equivalent under the action of Euclidean group. Counter-examples with 
closed curves and a discussion of the extend to which signature determines a smooth curve up to a 
group transformation can be found in |19|. 

4.3 Separating sets of invariants for afRne and projective actions 

In this section, we construct a differentially separating set of rational invariants for affine and 
projective actions. We will build them from classical invariants from differential geometry. We 



start with Euclidean curvature 
of the lowest order. Higher order Euclidean differential invariants are o 



a differential invariant 
Stained by differentiating 



^If X has a rational parameterization 7: R — > R^, then K\x and T\x are rational functions from R to R and, 
therefore, Sx is a rational map from R to R^. 

®For a parametric curve the formula k — ^^i^^^^iin can be used. 

'^The sign of k changes when a curve is reflected, rotated by vr radians or traced in the opposite direction. A 
rational function is invariant under the full Euclidean group. 
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the curvature with respect to the Euchdean arclength ds = \/l + [y^^^]'^ dx, i. e. 



d K 



d Ks 

ds ' ' 



d K 
d s 



^l + [y(l)]2dX^ ^' 

Affine and projective curvatures and infinitesimal arclengths are well known, and can be ex- 
pressed in terms of Euclidean invariants [8, 18 . In particular, 5^-curvature fi and infinitesimal 
5^-arclength da are expressed in terms of their Euclidean counterparts as follows: 



3 K {Kss + 3 K"^) 5 1/3 I 

a = ^ , da = K ' ds. 



(19) 



By considering the effects of scalings and reflections on 5^(2)-invariants, we obtain two lowest 
order ^(2)-invariants: 

Ka=^-^, T4=^. (20) 

They are of order 5 and 6, respectively, and are rational functions in jet variables]^ 

7-'^i2(3)-curvature ry and infinitesimal arclength dp are expressed in terms of their 5^-counterparts: 



V 



8/3 



dp = pl/^da. 



(21) 



The two lowest order rational P^£(3)-invariants are of differential order 7 and 8, respectively]^ 



Kp = 7] , Tp = rjp. 



(22) 



Theorem 24. ("differentially separating sets of affine and projective invariants. j 
According to Definitional^- 

1. The set Ij( = {K_a,TX\ given by (20) is differentially separating for the A{2)-action on 

2. The setX-pgc = {K-p^T-p} given by (22) is differentially separating for the VQC{?>)- action 
on M^. 

Proposition 25. I j[- exceptional algebraic curves are lines and parabolas. Xpgc- exceptional alge- 
braic curves are lines and conies. 



The proofs of Theorem |24| and Proposition 25 are given in Appendices 8.3 and 8.4, respectively. 

Corollary 26. An I j[- exceptional algebraic curve is not A(2) equivalent to a non I j[- exceptional 
algebraic curve. An Xpgc- exceptional algebraic curve is not VQC{?>) equivalent to a non Xpgc- 
exceptional algebraic curve. 



Theorem 24, in combination with Theorem 22 leads to a solution for the projective and the 
affine equivalence problems for curves. 



5 Algorithms and Examples 

In this section, we outline the algorithms for solving projection problems based on a combination of 
the projection criteria of Section [3] and the group equivalence criterion of Section [4} For simplicity, 
we treat the case of rational algebraic curves, while the algorithms covering arbitrary algebraic 



curves are presented in Appendix 8.5 



Explicit expressions of these invariants in terms of jet coordinates are given by formulas 
'Explicit expressions of these invariants in terms of jet coordinates are given by formulas 



(261 
(57i 



in Appendix 
in Appendix 
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5.1 Central projections. 

The following algorithm is based on the central projection criteria stated in Theorem [8j 
Algorithm 27. (Central projections.) 

INPUT: Rational maps T{s) = [zi{s), Z2{s), Z3{s)) , s G M andj{t) = {x{t), y{t)), t G R parametriz- 
ing algebraic curves ^ C M'^ and C M^. 

OUTPUT: The truth of the statement: 3[P] G CV, such that X = P{Z). 
STEPS: 



1. if X is a line or a conic then follow a special procedure (Algorithm 38 in Appendix 8.6), else 

2. evaluate VGC{3) -invariants (22) on X. Obtain two rational functions K-p\x andT-p\x oft; 

3. for arbitrary ci,C2,C3 G M define a curve parametrized by eJs) = ( ^^t^vt^, ^^r^YJ^h' 

4- evaluate VQC{'i) -invariants {22) on Zc- Obtain two rational functions K-p]-^ andT-pl^ of c 
and s; 

5. determine the truth of the statement: 

3c gM^ VsGM 3tGM K-pl^^ic, s) = K-plxit) and T-pl^^ic, s) = T-plxit). 

If the output is TRUE then, in many cases, we can, in addition to establishing the existence of 
ci , C2 , C3 in Step [5] of the algorithm, find at least one of such triplets explicitly. We then know that 
Z can be projected to A" by a projection centered at (— ci, — C2, — cs). We can also, in many cases, 
determine explicitly a transformation [A] G VGC{3) that maps X to Zc- We then know that Z can 
be projected to X by the projection [P] = [A][P^][B], where P^ and B are defined by Q. 



Example 28. We would like to decide if the spatial curve Z parametrized by T(s) = (s^, s^, s), 
s G M projects to any of the three given planar curves Xi, X2, X^ and A4 parametrized, respectively, 
by: 



t^ \ , . f t t^ 



For ci , C2 , C3 G M define a curve Zq parametrized by 6^(5) - — 1 ^ 1 5 

Note that Xi is parabola and so is 'PQC{2>)- exceptional. It is known that all planar conies are 
VQC{?>)- equivalent and so, from Theorem^ we know that Z can be projected on Xi if there exist 
c G M^, such that the curve Zc is a conic. This is obviously true for ci = C2 = C3 = 0. Indeed, one 
can check that Z can be projected to Xi by projection x = ^ , y = 

The curve X2 is not VQ C{3) -exceptional. It has a constant signature map: 



Following Algorithm 27, we need to decide whether there exist ci, 02,03 G M, such that K-pl^ (s) = 
T28(X)' ^''^^ '^'P\z = 0, Vs G M. This is, indeed, true for ci = C2 = and C3 = 1. We conclude 
that Z can be projected to X2. It is not difficult to find a possible projection: x = -^j^ , y = j^+j- 
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The curve has a non-constant signature map: 



9261 f -t^ + t 21 (t3 + 1)4 

Kv\x-,it) = — - —s- and T-plx-At) ~ 



50 (^3-1)8 ''-3v/ 10(t3-l)4- 

Evaluation of step^ of Algorithm 21 yields TRUE (one can, in fact, check that for ci = 0,C2 = 
0, C3 = 1 it is true that Kp\£ (s) = Kplx^is) and Tp\£ (s) = Kp\x2{s), Vs € Mj. We conclude that 
Z can be projected to X^. It is not difficult to determine a possible projection: x = -^j^ , y = jf^- 
It is important to observe that although Z can be projected to each of the planar curves Xi, X2, 
and X^, these planar curves are not VGC{3) -equivalent. This underscores an observation made in 
Remark^ 

The curve X4 has a constant signature map: K'p\x^{t) = andT'p\x^{t) = 0. Follow- 



ing Algorithm 21, we need to decide whether there exist ci, 02,03 G M, such that K-p\£ (s) 



and T'p\£^{s) = 0, Vs G M. Substitution of several values of s, yields a system of polynomial 
equations for ci, 02,03 G M that has no solutions. We conclude that there is no central projection 
from Z to X4. 

5.2 Parallel projections 



The following algorithm is based on the reduced parallel projection criteria stated in Corollary 10 



Algorithm 29. (Parallel projections.) 

INPUT: Rational maps T{s) = (zi(s), Z2(s), 2^3(5)), s G M and^{t) = (x(t), y{t)), t G '^parametriz- 
ing algebraic curves Z d'M^ and C M^. 

OUTPUT: The truth of the statement: 3[P] G VV, such that X = P{Z). 
STEPS: 



1. if X is a line or a parabola then follow a special procedure (Algorithm 40 in Appendix 8.6), 
else 



2. evaluate A{2) -invariants (20) on X. Obtain two rational functions Kj^\x and Kj^\x oft; 



3. define a curve Z parametrized by a(s) = {z2{s), 23(5)); 



4. evaluate A{2) -invariants (20) on Z. Obtain two rational functions Ka\z o^ndTj^l^ of s; 

5. determine the truth of the statement: 

VsGM 3teR Ka\z{s) = K^lxit) andT^lzis) =TA\x{t). 
If TRUE exit the procedure, else 

6. for arbitrary 6 G M, define a curve Z\, parametrized by /3b(s) = (zi(s) + 62:2(5), 23(5)); 



7. evaluate A{2) -invariants {20} on Zf). Obtain two rational functions Kj^l^ {b,s) andTj^l^ {b,s) 



of b and s; 

8. determine the truth of the statement: 

36 gM VsGM 3tGM Kj(\£^{b, s) = Kj^\x{t) and Tj^\£^{b, s) = Tj(\x{t). 

If TRUE exit the procedure, else 

9. for arbitrary a G M^, define a curve Za parametrized by 5a{s) = (21(5) + ai 23(5), 22 + «2 23(5)); 



14 



10. evaluate A{2) -invariants {20) on Za- Obtain two rational functions Kj[\£ (a, s) andTj[\£ {a,s) 
of ai, 02 and s; 

11. determine the truth of the statement: 

3a G Vs G M t G M K^\^^{a, s) = K_A\x{t) and Tji\^^{a, s) = Tji\x{t). 



Example 30. To decide whether the spatial curve Z parametrized by T{s) = Z2{s), z-^is)) = 

(s'* + 1, s^, s) , s G M, can be projected to X parametrized by ^{t) = (t , + t^) , t G M 6|/ a parallel 
projection, we start by noticing that X is not an A{2)- exceptional curve. Its signature map is given 
by non-constant rational functions: Tj[\x{t) = J^'^ap ^ i Kj[\x{t) = (i^'^^Q*^"^^^* +^) _ 



Following Algorithm 29. we first check whether X is A{2)- equivalent to Z parametrized by 
a{s) = (z2{s), Z3{s)) = {^P~s). The answer is no, since Z is an A{2)- exceptional curve (parabola) 
and X is not A{2) -exceptional. We next check whether there exists 6 G M such that X is .4(2)- 
equivalent to Zf, parametrized by f3b{s) = (^;i(s) + bz2{s), zs{s)) = (s^ + 1 + s). We evaluate 



invariants (20) on Zb: Tyi\^^{b,s) = ^^^^Sjj^jp^ , Kji^\^^{b,s) = ^ ^^'^(°li/j2)2' ^ ■ Since 
Ty(\2^{l, ,s) = Kj\\x and Tji\£^{l., s) = Tji\x, we conclude that Zb=i is A{2) -equivalent to X and, 
therefore, Z projects to X by a parallel projection. 

Example 31. We would like to decide if the spatial curve Z parametrized by (zi(s), Z2{s),z^{s)) = 
(s^ + s, s'^ — 3s^, s^) , s G M, projects to any of the three given planar curves Xi, X2 and X^ 
parametrized, respectively, by: 



7i(t) = (t^ + t, t^), j2{t) = {t^-t, t^ + t^), 73(0 =(- 



t^ 



1 + ^3)' (1 + ^3) 



Following Algorithm 29, we establish that Z parametrized by a{s) = (z2{s), ^^3(5)) and Z^ parametrized 
by Pb{s) = ( zi{s) + bz2{s), 23(5)), for all 6 G M, are not A{2) -equivalent to either of X 's. We then 
establish that Za parametrized by 6a{s) = {zi{s) + ai zs{s), Z2{s) + 02 -^3(5)) is A{2) -equivalent to 
Xi when ai = and 02 = 1/2 and is A{2) -equivalent to X2 when oi = and 02 = 0, but there are 
no real values of ai and 02 such that Za and X3 are A{2) -equivalent. 

We conclude that there are parallel projections of Z to both Xi and X2, but not to X3. Note that 
Xi and X2 are not A{2)- equivalent because (their affine signatures have different implicit equations). 
This underscores an observation made in Remark^ 



6 Projection of finite ordered sets (lists) of points 

In [l|[2], the authors present a solution to the problem of deciding whether or not there exists a 
parallel projection of a list Z = (z^, . . . , z"^) of m points in M'^ to a list X = (x^, . . . , x"^) of m 
points in M^, without finding a projection explicitly. They identify the lists Z and X with the 
elements of certain Grassmanian spaces and use Pliiker embedding of Grassmanians into projective 
spaces to explicitly define the algebraic variety that characterizes pairs of sets related by an prallel 
projection. 

We indicate here how our approach leads to an alternative solution for the projection problem 
for lists of points. Details of this adaptation appear in the dissertation 
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Figure 1: Projection problem for curves vs. projection problems for lists of points 



Theorem 32. (^central projection criteria for lists of points.) A given list Z = 
{z^, . . . , z™') of m points in with coordinates = {zi,Z2,z\), I = I . . .m, projects onto a given 
list X = (x^,...,x™) of m points in with coordinates x' = {x\y^), I = 1 . . .m, by a central 
projection if and only if there exist ci, C2, C3 G M and [A\ G VQC{'i), such that 

[x\ y\ 1]T = {A\{z{ + ci,zi + C2, 4 + C3f forl = l...m. (23) 



Theorem 33. (^parallel projection criteria for lists of points.) A given list Z = 
[z^, . . . , z™) of m points in with coordinates z' = {z[, zi^, ^3), I = 1 . . . m, projects onto a given 
list X = (x^, . . . , x™") ofm points in with coordinates x' = (a;', y'), I = 1 . . .m, by a parallel pro- 
jection if and only if there exist ci, C2 € an ordered triplet {i,j, k) G {(1, 2, 3), (1, 3, 2), (2, 3, 1)} 
and [A] G ^(2), such that 

[x\y\ if = [A] \z\ + ci zi, z\ + C2zllV forl = l... m. (24) 



The proofs of Theorems [32] and 33 are straightforward adaptations of the proofs of Theorems [8] 



and[9| The reduced parallel projection criteria for curves, given in Corollary |10[ is adapted to the 
finite lists in an analogous way. 

The central and the parallel projection problems for lists of m points is therefore reduced to 
a modification of the problems of equivalence of two lists of m points in PM^ under the action of 
VQC{2>) and ^(2) groups, respectively. A separating set of invariants for lists of m points in PM^ 
under ^(2)-action consists of ratios of certain areas and is listed, for instance, in Theorem 3.5 of 
20 . Similarly, a separating set of invariants for lists of m ordered points in PM^ under VQC{^)- 



action consists of cross-ratios of certain areas and is listed, for instance, in Theorem 3.10 in |20]. 
In the case of central projections we, therefore, obtain a system of polynomial equations on ci,C2 
and C3 that have solutions if and only if the given set Z projects to the given set X and analog of 
Algorithms |27| follows. The parallel projections are treated in a similar way. 

Figure [1] illustrates that a solution of the projection problem for lists of points does not provide 
an immediate solution to the discretization of the projection problem for curves. Indeed, let Z = 
(z^, . . . , z"^) be a discrete sampling of a spatial curve Z and X = (x^, . . . , x"^) be a discrete sampling 
of a planar curve X. It might be impossible to project the list Z onto X, even when the curve Z can 
be projected to the curve X . Some approaches to the discretization of the projection algorithms 
for curves are discussed in the next section. 

7 Directions of further research 

The projection criteria developed in Section |3] reduce the problem of object-image correspondence 
for curves under a projection from to to a variation of the group-equivalence problem for 
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curves in M^. We use differential signature construction [s] to address the group-equivalence prob- 
lem. In practical applications, curves are often given by samples of points. In this case, invari- 
ant numerical approximations of differential invariants presented in [3|[5] may be used to obtain 
signatures. Differential invariants and their approximations are highly sensitive to image pertur- 
bations and, therefore, are not practical in many situations. Other types of invariants, such as 



semi-differential (or joint) invariants 20,23 , integral invariants 10,13,21 and moment invariants 



17 are less sensitive to image perturbations and may be employed to solve the group-equivalence 



problem. 

One of the essential contributions of [l|[2j is the definition of an object/image distance between 
ordered sets of m points in and M^, such that the distance is zero if and only if these sets are 
related by a projection. Since, in practice, we are given only an approximate position of points, a 
"good" object/image distance provides a tool for deciding whether a given set of points in is 
a good approximation of a projection of a given set of points in M^. Defining such object/image 
distance in the case of curves is an important direction of further research. 



8 Appendix 

8.1 Explicit formulas for invariants 

Let 

Ai = 3y(^) _ 5 [y(3)]2 ^^^^ ^ 9^(5) [y(2)]2 _ 45^(4) ^(3) ^(2) ^ 40 [^(3)]= 



The explicit formulas for the differentially separating set of rational ^(2)-invariants ( 20 ) are: 



Ka 



Ta 



(Ai)3' 

9 ^(6) [^(2)]3 _ 63 ^(5) ^(3) [^(2)j2 _ 45 [^(4)]2 [^(2)p ^ 255 y ^ [y(3)]^ y(^) - 160 [j/(3)]^ 

(Ai)2 



(25) 



(26) 



The explicit formulas for the differentially separating set of rational T'^£(3)-invariants (22) are: 
729 



8(A2 



18y(^) [y(2)]4 A2-I89 



(6)12 



2)16 



(27) 



+ 126 (9 y(5) y{3) y(2) ^ ^5 [y{4)]2 ^(2) _ 35 y(4) [y{3)]2^ 



189 (4[y(^)]^ + 15y(2)y(^)) 



+ 210y(^)y(3)[y(2)]2 



(63 [y lyi-)^ _ 60y(4) [y(3)]2 ^(2) + 32 [^(3)]. 



- 525 y (9 [y^'^f [y^'^f + 15 [y^'^f [y^'^f [y^'^]' - 60 y(^) [y^'^]' y^'^ + 64 [y('^f) 
+ 11200 [y(3)]^ 
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2(A2 



- 8y(7) A2 (9y(6) [y{2)]3 _ 36^(5) ^(3) [y(2)]2 _ 45 [y(4)p [y{2)]2 ^ ^20y(4) [y(3)]2y{2) _ 40 [y(3)]^ 

+ 504 [y(6)]3 [y(2)]5 _ 504 [y(6)]2 [y(2)]3 (9 ^(5) y(3) y(2) ^ 15 [y(4)]2 ^(2) _ 35 ^(4) [y(3)]2) 

28 (432 [y(5)]2 [y(3)]2 [y(2)]3 + 243 \y^'>\^ y ^ [y(2)]4 _ igoo [y(3)]3 [y{2)]2 

240y(5) [y{3)]5 y{2) ^ 54Qy(5) [y{4)]2 [y(3)] [y(2)]3 ^ ggQQ [y(4)]2 [^(3)]4 ^(2) _ 2000y(4) [y{3)]6 
5175 [y(3)]2 [y(2)]2 + 1350 [y(4)]4 [y(2)]3) _ 3335 {y^^\' [y(2)]4 

252 [y(5)]3y(3) [y(2)]2 (9^(4) ^(2) _ ^gg [^(3)^2) _ 35340 [y(5)]2 [y(3)]6 

- 630 [y [^(2)] (59 [y(4)]2 [^(2)^2 _ ^gg [y(3)]4 _ 153 y(4) [y(3)]2 [^(2)]) 

+ 2100 y(5) [y(4)]2 y^^) (72 [y(3)]^ + 63 [y^^)]^ [y(2)]2 _ 193 y^ [y(3)]2 ^(2)) 

- 7875 [yW]^ (8 {y^'^f [y(2)]2 - 22y(^) \y^^f \y^^\ + 9 [y^^)]^] 



+ 



+ 



Remark 34. // Ai is zero at more than at a finite number of points of an algebraic curve X then 
X is either a line or a parabola, and then Ai = 0. If ^2 is zero in more than at a finite number of 
points of X then X is either a line or a conic and then A2 = 0. 

8.2 Proof of Theorem [22l 

Direction =^ follows immediately from the definition of invariants. Below we prove <^^. We notice 
that there are two cases. Either K\x^ and K\x2 are constant maps on Xi and ^2, respectively, and 
these maps take the same value. Otherwise both K\xj^ and K\x2 are non-constant rational maps 
on Xi and X2, respectively. 

Case 1: There exists c E M such that K\xi (pi) = c and Kx2 (P2) = c for all pi G Xi and for all P2 G ^^2- 
Since Xi and X2 are non-exceptional, we may fix Xc-regular points pi = (xi,yi) G Xi and 
P2 = {x2,y2) € X2- Then, due to separation property of the invariant K , 3g & G such that 



j^^"^(pi) = g ■ [j^2^(P2)]- consider a new algebraic curve X3 = g ■ X2. Then due to (16), 
we have 

fx-\pi) = fx:\pi) =■■ P^'~'^ ■ (28) 



Since pi is a X-regular point of Xi, it follows from ( |28| ) that it is also a X-regular point 
of X3 and, in particular, is non-singular. Let Fi{x,y) = and i<3(a;,y) = be implicit 
equations of Xi and X3 repsectively. We may assume that ^^(pi) 7^ and ^^(pi) / 
(otherwise, ^?^(pi) 7^ and ^^^(pi) 7^ and we may use a similar argument). Then, there 
exist functions fi{x) and /3(x), analytic on an interval / B xi, such that Fi{x, fi{x)) = 
and F3{x, f3{x)) = for x G /i. 

Functions y = fi{x) and y = fsix) are local analytic solutions of differential equation 

K(x,y,yW,...,y(^-i)) = c (29) 



with the same initial condition f[''\xi) = fl^''\xi), A; = 0, . . . , r — 1 prescribed by (28). From 



the X-regularity of pi, we have that Qy(^-i) / so (29) can be solved for y^*" 



1) 



y(^-i)=i/(x,y,y«,...,y(^-2)), (30) 
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where function H is smooth in a neighborhood p^''"^) G J^~^. From the uniqueness theorem 
for the solutions of ODEs, it foUows that fi{x) = /3(x) on an interval I 3 xi. Since Xi and 
A's are irreducible algebraic curves it follows that Xi = X3. Therefore, Xi = g ■ X2- 

Case 2: K\x-^ and K\x2 are non-constant rational maps. Then S^i = is a one-dimesional set that 
we will denote S. Let S{x,t) = be the implicit equation for S (see Remark 20). We know 
that ^{k,t) 7^ for all but finite number of values (x, r), because, otherwise, K\x-^ and 
K\x2 are constant maps and the set S is zero-dimenional. Therefore, since the curves are 
non-exceptional, there exists X-regular points pi = (xi,yi) G Xi and p2 = (2:2,2/2) £ X2 such 
that 



^Ui(pi) = K\x2{V2) =■■ xo, T\x^ipi) = T|a'2(p2) 



dS 

--: To and ^(xo,ro) / 0. (31) 



Due to separation property of the set Iq = {K,T}, 3g £ G such that j^j(pi) = g ■ [j^2(P2)]- 
We consider a new algebraic curve X^ = g ■ X2. Then due to (16), we have 



(r) 



From (31), (32) and X-regularity of the point pi G X\ it follows that 

dT 



xo, T(p 



tq and 



^0 



(32) 



(33) 



Since pi is a X- regular point of Xi, it follows from (32) that it is also a X-regular point of X-- 



and, in particular, is non-singular. Let Fi{x, y) = and Fs{x, y) = he implicit equations of 
Xi and X3 repsectively. We may assume that 7^ and 7^ (otherwise, ^ and 

7^ and we may use a similar argument). Then, there exist functions /i(x) and fs^x), 
analytic on an interval I B xi, such that Fi{x, fi{x)) = and F3{x, fsix)) = for x G Ii. 

Then functions y = fi{x) and y = fz{x) are local analytic solutions of differential equation 



S (k{x, y, , . . . , y(^-i) ) , T{x, y, y^D , . . . , yM)) = 



(34) 



with the same initial condition f'^\xi) = flf^'{xi), k = 0, . . . ,r, dictated by (32). 



,ir). 



y 



(r) 



Since ^(xo,ro) 7^ and -7^^ ^ ^ 7^ (see (31) and (33)), equation (34) can be solved for 



H{x,y,y 



(1) 



<y 



(r-l)^ 



(35) 



where function H is smooth in a neighborhood p^''^ G J*". From the uniqueness theorem for 
the solutions of ODE it follows that fi{x) = fsix) on an interval I 3 xi. Since Xi and X3 
are irreducible algebraic curves it follows that Xi = X3. Therefore, Xi = g ■ X2. 



□ 
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8.3 Proof of Theorem I 

1. We note that dim ^(2) = 6. We will prove the separation property of X4 = {Ky\^,Ty\} given 



by ( 26 ) on a Zariski open subset 

= {p(6) G J6 I ^ and Sy^^) yW _ 5 [yi3)f q} 

of and the separation property of Ky^ on = 7t^{W^). An affine transformation can be 
written as a product of a Euclidean, an upper triangular area preserving linear transformation, 
a scaling and a reflection. 



A 




f 

1 

1 




where + S'^ = 1, e = ±1 and h ^ 0. 
A Euclidean transformation with 

1 y^^^ y^^'^ y + X y^^'^ x 



Vi+WW VT+WW VT+WW yr+F^ 



(36) 



brings a point p(^) = (x, y, y^-*^), y(^), y(^), y^^)) to a point p^^"* = {0,0,0,y^\y^i^\y[^\y[^^), 



where 7^ if and only if y^'^' ^ 0. An upper triangular transformation with e = 
[yf and / = Jajis/s brings the point p^^^ to a point p^^^* = (0, 0, 0, 1, 0, y2\yf^), where 



7^ if and only if Sy'-'^-' y*-^^ — 5 [y^'^^]^ 7^ 0. A scaling transformation with h = [y2^^]3 
brings pf ^ to a point p^^^ = (0, 0, 0, 1, 0, 1, yf''). Finally, if yf^ < 0, we can apply a reflection 
(e = -1) to bring p^^^ to pf ^ = (0,0,0, 1,0, l,yf with y^ > 0. 

By the above argument, for any p(^\ q^^^ G W^, there exist A, B £ -4.(2) such that p^^^ = 

A ■ p(5) = (0,0,0,l,0,l,yP) and qf^ = B ■ q^^) = (0, 0, 0, 1, 0, 1, y^^^), where yf^ > and 

(5) 

y2 ^0. Then since Kj[ is invariant: 

Ka{p^'^) = Ka{p?) = 3 [yP]2 and K^ci^'^) = K^cL?) = 3 [y?^]^ 
Since yf'^ > and y^^ > 

Ka{p^'^) = KA^^'^) ^ y? = yf^ ^ = ^ q^^^ = BA-^ . p(^). 

and therefore is separating on W^. 

Similarly, for any p('^),q('^) G W^, there exist A, B £ A{2) such that pf^ = A ■ p(^) = 
(0, 0,0,1, 0,l,yP,yf) and qf^ = i? • q(6) = (0, 0, 0, 1, 0, 1, yf \ yf ), where yf^ > and 
2/2^^ > 0. Then since and are invariant: 

(6)n _ Qr„.(5)l2 r i^.r„(6)N _ _ o r„.(5)i2 



i^^(p(^)) = i^^(pr) = 3[yff i^^(q(6)) = = 3 [y 



^2 



T^(p(^)) = T^(pf^) = ""'1 r^(q(^)) = T^(qf)) = yf 
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Then 



i^^(p(')) = A'^(q(^)) { yf = ^ (6) _ (6)^^(6) _ 



T^(p(^)) = T^(q(^)) = ^Pi -qi -p 

and thus the separation property of {Ky\^, T4} on is proved. 

2. We note that dimVQC^S) = 8. We wih prove the separation property of 2p = {K-pjT-p} 
given by (27) on a Zariski open subset 

= {p(8) e J« I ^ and Oy^^) [y(2)]2 _ 45y(4) ^(3) ^(2) ^ [^(3)^3 o| 
of and the separation property of Kp on W''' = '^^{W^). 

Observe that any element in ^ G VQC{3) can be written as the following product: 

/ 



A 




1 




where + = 1, e 7^ and g ^ 0. In the first part of the proof, we have shown 
that the two matrices on the right can bring a point p*^^) G to the point p^^'' = 
(0, 0, 0, 1, 0, y[^\ . . . , 2/1^^)- Observe that is invariant under these transformations. There- 
fore p^^^ G and hence y^^^ 7^ 0. The last transformation with 

, _ r,(5)ii/3 , _ 1 5[yf]^ , ^ _ 1 [y'^'?-i0y'^'[yi''?-nyP?y'^' + 25[y['^]^ 

g-[y^ I ,n- ana z - 

brings to = (0, 0, 0, 1, 0, 0, 1, 0, 7/2 1 ?/2 )• direct computation shows that Kp{p2 ) = 
[2/2'''^]^ and Tp(p2^^) = 3y^^ — We observe that y^^ and y^2^ can be uniquely determined 
from the values of invariants K-p and T-p and therefore we can complete the proof of the 
separation property of Kp on VF^ and the separation property of {Kp,Tp} on by an 
argument similar to the one presented in part 1 of the proof. 

□ 



8.4 Proof of Proposition | 

In the affine case, we note that j%{p) G \ and j%{p) G \ if an only if k\x{p) = 
or fi\xip) = 0. If k|;\:'(p) = for more than finite number of points on X, then is a line, if 
f^\x{p) — foi' more than finite number of points on X then it is a parabola (see |T2|). From 
the exphcit formulas (26) we see that / and |% / for aU p(5) G and all 

^ ayWJ p{5) ' dyy'^) p(6) ' ^ 

p(6) g . Therefore, if an algebraic curve is not a line or a parabola it is {/C^, 73i}-regular. In 
the projective case, we note that ij(p) S \ and j%{p) G \ W'^ if an only if k\x{p) = or 
^ia\x{V') = 0- If /^qIa'Cp) = for more than finite number of points on X then A:" is a conic (see |12|). 



From the exphcit formulas \27\, we see that 0^ ^ and / for all pC'') G VF''' and 



-,(8) 



all p^*^^ G ly^. Therefore, if an algebraic curve is not a line or a conic it is {/C-p, 7p}-regular. □ 



21 



8.5 Projection algorithm for arbitrary algebraic curves 



In this section we no longer assume that algebraic curves X and Z have a rational parameterization. 
We continue to assume that Z is not a straight line. 

Algorithm 35. (Central projections.) 

INPUT: An irreducible polynomial F(x, y) whose zero set is an algebraic curve X and a prime 
idea H C M[zi, Z2, -^s] defining algebraic curve Z <ZW^ . 



OUTPUT: The truth of the statement: 3[P] G CV, such that X = P{Z) 
STEPS: 



1. if X is a line or a conic then follow a special procedure (see Appendix 8.6), else 

2. evaluate VGC{3) -invariants (22) on X, using (15). Obtain two rational functions K'p\x and 
T-p\x of X and y; 



3. use an elimination algorithm (see Remark 36) to compute an irreducible polynomial 
G G M[ci, C2, C3, x, y], which for a fixed c G defines the curve 

p2 . 



Z3+C3 ' 



Z2+C2 



(^1,^2,^3) & Z} C 



4. evaluate VQC{3) -invariants (22) on Zc, using (15). Obtain two rational functions K-plf 
Tp\£^ of c, X andy; 

5. determine the truth of the statement: 

3c G y{x,y) e X 3(x,y)GM^ G{c,x,y) = 0, K-p\£^{c, x,y) = K-p\xix,y) and 

T'p\zic, x,y) = Tp\x{x,y). 



and 



Remark 36. To execute Step^ of Algorithm 35, we define an ideal r. = H -\- (^x (2:3 + C3) — (zi + 
ci)) y (-2^3 + C3) — (-22 + C2), 5 (-^3 + C3) — 1) C M[c, X, y, zi, Z2, Z3, S\ . An elimination algorithm Can 
be used to compute the intersection H = H n M[c, The zero set of H is the closure of the 

image of the rational map (p, defined by (/)(ci, C2, C3, zi, Z2, ^3) = (ci, 02,03, f|$^) > from the 

irreducible variety cp: x Z to M^, and, therefore, is an irreducible variety. Since the dimension 
of this variety equals to 4, it is defined by a single irreducible polynomial G{c,x,y). 

Algorithm 37. (Parallel projections.) 

INPUT: An irreducible polynomial F{x,y) whose zero set is an algebraic curve X and a prime 
idealH C M[zi, 22, -23] defining algebraic curve ^ C M^. 

OUTPUT: The truth of the statement: 3[P] G VV, such that X = P{Z). 

STEPS: 



1. if X is a line or a parabola then follow a special procedure (see Appendix 8.6), else 

2. evaluate A{2) -invariants (20) on X , using (15). Obtain two rational functions Kyx\x o-nd 
Ta\x ofx and y; 



^"^Generically, H has two generators, although twisted cubic provides a famous counter-example. The question 
whether or not any algebraic curve in (or C'^) is a set-theoretic intersection of two algebraic surfaces remains open 
151. 
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3. use an elimination algorithm to compute an irreducible polynomial G G M[x,y], which defines 
the curve Z = ^{z2, 23) (2:1,^2,-23) G C M^; 

4. evaluate A{2) -invariants (20) on Z using (15). Obtain two rational functions K_a\£ and Tj^l^ 
of X and y; 

5. determine the truth of the statement: 

y{x,y)eX G = 0, Kji(\^{x,y) = Kj\\x{x,y) and Tji\2{x,y) = Tj^\x{x,y). 

If TRUE exit the procedure, else 

6. use an elimination algorithm to compute an irreducible polynomial G G M[6, which for 
a fixed 6 G M defines the curve Z}y = |( 21 + 62:2; ^3) 



(^1,22,^3) G C 



7. evaluate A{2) -invariants (20) on Z^, using (15). Obtain two rational functions KaIz,, '^^'^ 



TaIz^ 0/6, X and y; 

8. determine the truth of the statement: 

36 gM V(x,y)GA' 3(f,y)GM2 G(6,x,y) = 0, KA\z^{b, x,y) = Kj^\x{x,y) and 

If TRUE exit the procedure, else 

9. use an elimination algorithm to compute an irreducible polynomial G G M[ai, 02, x, y], which 

{zi,Z2,Zs) £ Z^ CM2; 



for a fixed a G defines the curve = |( zi + ai Z3, Z2 + 02 -23) 



10. evaluate A{2) -invariants (20) on Za using (15). Obtain two rational functions Kj^\£ and 
TAlza '^2; X andy; 

11. determine the truth of the statement: 

3a G y{x,y) £ X 3(£,y)GM^ G{ai,a2,x,y) = 0, Kyi\^^{ai,a2, x,y) = Kj[\x{x,y) and 

^^Ij, (01,02, x,y) = T_A\x{x,y). 

8.6 Algorithms for exceptional curves 

We assume here that Z is a rational curve. The case when Z is not rational can be treated similarly. 
Algorithm 38. (Central projections to lines and conics.) 

INPUT: Rational maps r(s) = (2:1(5), 2:2(5), 2:3(5)), s G M and^{t) = (x(t), y{t)), t G M parametriz- 
ing algebraic curves 2^ C and X C M?, such that A2\x{t) = (see (25)). 

OUTPUT: The truth of the statement: 3[P] G CV, such that X = P{Z). 

STEPS: 

1. ifxy — yx = and (F x F) • r = return TRUE and exit; 

2. ifxy — yx = and (F x F) • F 7^ return FALSE and exit; 

3. for arbitrary ci, 02,03 G M define a curve Zc parametrized by ec(s) = , 
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4- evaluate A2IJ using (25). Obtain a rational function A2I2 of c and s. 
5. determine the truth of the statement: 3c G Vs € M A2IJ (c, s) = 0. 

Remark 39. The first two steps of the above algorithm rely on the fact that a space curve can be 
projected to a line if and only if the space curve is coplanar. 

Algorithm 40. (Parallel projections to lines and parabolas.) 

INPUT: Rational maps T{s) = (zi(s), Z2(s), 2:3(5)), s E M and^{t) = {x{t), y{t)), t G R parametriz- 
ing algebraic curves Z cM.^ and X C M'^, such that /S.i\x{t) = (see (25)). 

OUTPUT: The truth of the statement: 3[P] G CV, such that X = P{Z). 

STEPS: 

1. ifxy — yx = and (F x F) • r = return TRUE and exit; 

2. ifxy — yx = and (F x F) • F 7^ return FALSE and exit; 

3. define a curve Z parametrized by a(s) = (2:2(5), 23(5)); 

4-. evaluate Ail^ using (25). Obtain a rational function Ai|^ of s; 

5. determine the truth of the statement: Vs G M Ai|j(s) = 0. 
// TRUE exit the procedure, else 

6. for arbitrary 6 G M, define a curve Z\y parametrized by /3f,(s) = (21(5) + 622(5), 23(5)); 

7. evaluate Ail^^^ using (25). Obtain a rational function Ai|j^ of b and s; 

8. determine the truth of the statement: 36 G M Vs G M Ai|j^(6, s) = 0. 
// TRUE exit the procedure, else 

9. for arbitrary a G M^, define a curve Za parametrized by 5a{s) = (21(5) + ai 23(5), 22 + 02 23(5)); 
10. evaluate Ai\£ using (25). Obtain a rational function Ai\£ ofai,a2 and s; 



11. determine the truth of the statement: 3a G Vs G M Ai|j (a, s) = 0. 
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